Elliptic flow of thermal photons and formation time of quark gluon plasma at RHIC 
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We calculate the elliptic flow of thermal photons from Au+Au collisions at RHIC energies for a 
range of values for the formation time to but a fixed entropy (or particle rapidity density). The 
results are found to be quite sensitive to tq. The V2 for photons decreases as to decreases and admits 
a larger contribution from the QGP phase which has a smaller V2- The elliptic flow coefficient for 
hadrons, on the other hand, is only marginally dependent on tq. 
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It is by now generally believed that we are witness- 
ing the birth of quark gluon plasma in relativistic colli- 
sion of gold nuclei at RHIC energies. This view is sup- 
ported by the clear display of jet-quenching elliptic 
flow of particles in non-central collisions [3, Ij] heralding 
the on-set of collectivity and conversion of initial spa- 
tial anisotropy to azimuthal anisotropy of their momenta, 
and radiation of thermal photons [5| . Now the next set 
of questions occupy the central stage in this endeavour. 
What is the temperature of the plasma? Is it in ther- 
mal and chemical equilibrium? What is the formation 
time To, beyond which we could possibly use the power- 
ful methods of hydrodynamics to describe the evolution 
of the system? 

The question of tq is also related to the process of ther- 
malization. In the simplest treatment one assumes that 
the partons produced in the collision have an average 
energy (E) and thus tq ~ 1/ (E) . It is a different ques- 
tion whether these are thermalized 6] . One may use one 
of the many mini-jet based models and assume that the 
energy and momentum of the partons are "some- how" 
redistributed and the plasma is formed at the deposited 
energy density [3] . A further refinement may include use 
of a parton cascade model [1, 0, [l3| , a self-screened par- 
ton cascade model [lH , or parton saturation models 
to get a time tq ~ l/pr or as/ Qs (Ref. at which one 
assumes the plasma to be thermally equilibrated. The 
chemical equilibration can also be modeled using master 
equations 0, 0] . Even though discussed in literature for 
quite some time [lH, the role of plasma instabilities in 
thermalizing the plasma is being explored in detail only 
now (see Ref. 16]). 

In view of the several treatments and perhaps compli- 
mentary models available in the literature, a more direct 
measurement of the formation time tq would be very de- 
sirable. In the present work we argue that a experimen- 
tal determination of the V2 of thermal photons for non- 
central collisions can help us determine this value rather 
accurately. 

Elliptic flow of thermal photons and dileptons 17 , 1^ 
has recently been proposed as a powerful tool for deter- 
mining the azimuthal anisotropy of momenta of partons 
soon after the formation of the quark- gluon plasma. The 
V2 for hadrons though seeded in the early QGP phase (as 
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FIG. 1: (Colour on-line) Spectrum (upper panel) and «2 
(lower panel) for primary pions from Au-|-Au collisions at 200 
GeV/A for 6 = 6 fm for different initial times. 



suggested by the re-combination models [19| and quark- 
number scaling of hadronic V2) determines this indirectly. 

We start with the assumption that a chemically and 
thermally equilibrated quark gluon plasma is produced 
at time tq beyond which we trace its evolution using an 
ideal hydrodynamics. We closely follow the procedure 
outlined earlier [l3,[l3| in setting up the initial conditions 
and assume that the entropy density s(ro, x, y, 6) is given 
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by: 



s{To,x,y,b) = K [anw{x,y,b) + (1 - a)nb{x,y,b)] (1) 

where k is a constant, is number of wounded- nucleons 
(participants) at impact parameter 6, and rib is the 
number of binary coUisions estimated using the Glauber 
modeL In addition, a — 0.75 is a constant which controls 
the relative contributions of soft and hard processes. We 
choose K, such that so — s{tq — 0.6 fm, x — 0, y = 
0, b ^ 0)= 117 fm^^ for the collision of gold nuclei 
at ^/s = 200 A GeV. We adjust k to get the starting 
values for sq at other values of tq such that sqTq (or the 
particle rapidity density [21| ) is a constant, which implies 
identical entropy for all the cases. Note that the profile 
of the entropy distribution is uniquely determined by the 
above equation. We employ an impact parameter de- 
pendent, boost invariant hydrodynamics [22|, which has 
been used extensively to explore hadron production and 
elliptic flow of hadrons as well as photons 17[ and dilep- 
tons [11. 
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FIG. 2: (Colour on-line) Same as Fig. [T]for rho-mesons. 

The production of thermal photons is calculated by 
folding the history of the evolution of the system with 
the rate for the production of photons from the quark 
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FIG. 3: (Colour on-line) Effect of changing the freeze-out 
energy density on the spectra (upper panel) and V2 (lower 
panel) for tt and p mesons. The plot for V2 for the p-mesons 
is shifted by 0.5 GeV/c along the x-axis. 



matter and the hadronic matter. We use the complete 
leading-order results for the production of photons from 
the QGP [2^ and the latest results for the radiation of 
photons from a hot hadronic gas [24] . As mentioned ear- 
lier, the equation of state (EOS Q 22]) incorporating a 
phase transition to quark gluon plasma at T w 164 MeV, 
and a resonance gas for the hadronic phase below the 
energy density of 0.45 GeV/fm~'^ is used to describe the 
evolution. The mixed phase is described using Maxwell's 
construction. The freeze-out is assumed to take place at 
e/ = 0.075 GeV/fm^. The spectrum of hadrons is ob- 
tained using the Cooper-Fry formulation. 

As a first step we show the calculated spectra of pri- 
mary pions and rho mesons along with the differential 
elliptic flow for a typical impact parameter of 6 fm for 
Au-|-Au collisions at RHIC energies. We see that if the 
initial time is large, the inverse slope of the pion spec- 
trum (Fig. [TJ upper panel) is smaller due to a late start 
of flow (all the results are for zero rapidity) . However the 
deviation in the spectrum up to about 1.5 GeV/c is not 
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substantial, even though the spectrum at pt ~ 4 GeV/c 
changes by about a factor of three. It is very interesting 
to note that the differential elliptic flow is almost inde- 
pendent of the initial time across the range of px values 
considered. Similar results are obtained for the primary 
rho mesons, even though the change in the spectra is 
more noticeable as the rho meson (see Fig. [2]) being heav- 
ier is more strongly affected by the flow. We have verified 
that similar results are obtained for all other particles, 
included in the (nearly complete) list of hadrons in our 
calculations. The similarity of spectra for the primaries 
ensures that the spectra of particles after the resonance 
decays are accounted for, will remain similar. 

One may argue that the freeze-out at a given (rather 
low) energy density represents an over-simplification 
of the decoupling of the hadrons from the interacting 
medium. Thus one could consider hadrons decoupling at 
a sufficiently large temperature or energy density from 
the hydrodynamic evolution and populating a hadronic 
cascade from which they undergo a freeze-out depending 
on their interaction cross-sections [i^, . While a full 
calculation along these lines would be welcome, we can 
get a feeling of the outcome by estimating the effect of 
the freeze-out density on the spectra and the V2 of the 
hadrons using hydrodynamics. 

In Fig. [3l we give our results for the pt distributions 
and differential elliptic flow, V2{pt), for pions and p- 
mesons for energy- densities of 0.045 GeV/fm^ and 0.135 
GeV/fm'^, which correspond to freeze-out temperatures 
of about 100 and 140 MeV respectively. We note that the 
Pt distribution continues to evolve as the system cools 
from 140 MeV to 100 MeV, due to the radial flow. How- 
ever, the elliptic flow is seen to have essentially acquired 
its final value when the temperature is still large. Thus 
we feel that if these hadrons are fed into a hadronic cas- 
cade the results for the V2 will remain unaltered. 

We conclude thus, that it would not be easy to deter- 
mine To on the basis of either spectra or V2 for hadrons. 
The calculated differences in the spectra at very large val- 
ues of Pt may not be quite useful, as the hydrodynamic 
description itself may not hold at ^ 2 GeV/c. 

Next we calculate the thermal photon production for 
the same system (Fig. |4]) using tq varying from 0.2 fm/c 
to 1 fm/c. It is seen that the contributions from the 
hadronic matter for different starting times are only 
marginally different at lower pr, though at 6 GeV/c it 
increases by an order of magnitude, as tq drops from 1.0 
fm/c to 0.2 fm/c, due to an increased flow for lower initial 
times. However, the smaller initial times also admit much 
larger initial temperatures and the radiation of high pT 
photons from the quark matter increases substantially. 
Thus, while the total production of thermal photons for 
Pt < 1 GeV/c is not strongly affected, aX pt — Q GeV/c, 
it increases by more than four orders of magnitude. This 
is known for quite some time in the literature. 

The differential elliptical flow for thermal photons 
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FIG. 4: (Colour on-line) Spectrum of photons radiated from 
hadronic matter alone (upper panel) and the sum (lower 
panel) for different initial times. 



shows a much more complex and rich structure (Fig. [5]). 
The V2 for photons from the quark matter is much smaller 
than that for photons from the hadronic matter [l^ . 
This, happens as the anisotropics in momenta of the 
particles producing photons take some time to develop 
and the photons from the quark matter measure these 
anisotropics for early times. Note also the peak around 
0.5 GeV/c due to the dominating role of the reaction 
TT/d — > 7r7 [l3| beyond about 0.5 GeV/c. We also see an 
interesting trend; the V2 for thermal photons from the 
quark matter increases with the increase in the initial 
time To, while those for the photons from the hadronic 
matter decrease with increase in the initial time tq. This, 
we feel, happens due to a competition between larger ini- 
tial (spatial) anisotropics and temperatures for smaller tq 
and a start of cooling and equalization of pressure gradi- 
ents at an earlier time. 

The results for the final V2 for thermal photons (solid 
curves. Fig. [5|) reveal a large sensitivity to the initial 
time To, for pT greater than about 1.5 GeV/c. This is 
brought about by the increasing contribution of photons 
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FIG. 5: (Colour on-line) The differential elliptic flow of ther- 
mal photons (solid curves) for different initial time tq. The 
results for emissions from the hadronic matter (long-dashed) 
and quark matter (short dashed) are also shown. 
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FIG. 6: (Colour on-line) Single photons for 10-20% most cen- 
tral colhsions for Au-|-Au system at 200 GeV/A. NLO pQCD 
results for prompt photons (see text) and thermal photons for 
different formation times are also shown. 



from quark-matter (with smaller 112) as the time tq de- 
creases. With increasing tq the relative contribution of 
thermal photons from the hadronic matter increases. As 
these have larger W2, the over-all vi becomes larger. 

Before closing we show our results for single photon 
production and compare it with the measurements by 
the PHENIX experiment d for the 10-20% most cen- 
tral collisions for Au-|-Au system at 200 GeV/A (Fig.^. 
The prompt photon production has been calculated using 
NLO pQCD 27[ and choosing the factorization, renor- 
malization, and fragmentation scales as equal to pt/2. 
The effects of (impact parameter dependent) shadow- 
ing (2^ and iso-spin is explicitly accounted for. CTEQ5M 
structure functions are used. The thermal photon pro- 
duction for different initial times is also shown. We see 



that the NLO pQCD calculations nearly fully account for 
the production of photons at pt > 3 GeV/c. 

We would like to add though that the fragmentation 
contribution to the prompt photons (about 30% for -pT > 
3 GeV) will be suppressed due to jet-quenching. Assum- 
ing a suppression by a factor of about 2 (for quark jets 
which fragment into photons) this would mean a reduc- 
tion by « 15% in the pQCD yield compared to what is 
shown here. However, this would be somewhat off-set by 
the production of photons due to passage of jets through 
the QGP [29i]. A more detailed and complete calculation 
will consider all these sources [s^l- We postpone that to 
a future publication as the present emphasis is on the az- 
imuthal anisotropy of thermal photons, which is rather 
negligible for other direct photons. Similarly, one may 
consider the pre-equilibrium contribution as well |31| , if 
To is assumed to be large. These contributions will not 
show any V2 as they are not subjected to collectivity. 

If we are unable to separate these other sources of 
photons by, say, either by tagging or by calculations, 
then the V2 given here will have to be moderated as 
(vf X iVth + X iVnon-th)/(^th + N^on-tu), whcre 

'th' stands for 'thermal'. However as is quite 

small, the difference in the resulting V2 for single pho- 
tons will still persist, being large for large tq. In this 
connection, the success of the NLO pQCD in predicting 
the prompt contribution can be used with a great advan- 
tage as it can be subtracted from the direct photons to 
greatly increase the sensitivity of the remainder to tq . 

In order to get a more quantitative idea about 
these, we draw the attention of the readers to FIG. 6 
of Ref. [3^. First of all, as remarked earlier, the 
Compton+ annihilation part of the prompt photons does 
not contribute to azimuthal anisotropy. The prompt- 
fragmentation part, again as suggested earlier, will show 
a marginal positive V2, due to the reaction plane depen- 
dence of energy loss suffered by the out-going quarks be- 
fore fragmenting. This according to Ref. [sO] is less than 
about -fl%. The jet conversion photons give rise to very 
small (less than —1%) V2- The resulting V2 due to these 
two processes is essentially zero [sBj. This, along with 
the possibility of an accurate determination of prompt 
photons using NLO pQCD as seen here holds out the 
promise of determining the sensitivity of elliptic flow of 
thermal photons to tq. 

An experimental verification of these predictions will 
give a very useful information about the formation time 
and also validate our ideas about the evolution of the 
hot and dense system in relativistic heavy ion collisions. 
A theoretical reanalysis of single photon production for 
Pb+Pb collisions at the CERN SPS has revealed similar 
sensitivity to the formation time |32i | . 

In brief, we have calculated the thermal photon pro- 
duction for non-central collisions of gold nuclei at RHIC 
energies. The initial time is varied from 0.2 fm/c to 
1 fm/c, keeping the total entropy fixed and assuming 
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formation of quark gluon plasma in thermal and chemi- 
cal equilibrium, undergoing expansion and cooling. The 
V2{pt) for hadrons is seen to be only marginally affected 
by the variations in the formation time. The spectra 
as well as the V2{pt) of thermal photons are seen to be 
strongly affected by the formation time, as smaller ini- 
tial times admit larger initial temperatures and larger 
contributions from the QGP phase, when the momen- 
tum anisotropy is not fully developed. This could prove 
useful in determining the formation time tq. 

We thank Sangyong Jeon for his help in setting up 
the prompt photon calculation program of Patrick Au- 
renche's group. 
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